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Abstract 

An algorithm to characterize collective motion is presented, with the in- 
troduction of "collective Lyapunov exponent", as the orbital instability at 
a macroscopic level. By applying the algorithm to a globally coupled map, 
existence of low-dimensional collective chaos is confirmed, where the scale 
of (high-dimensional) microscopic chaos is separated from the macroscopic 

motion, and the scale approaches zero in the thermodynamic limit. 
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Low- dimensional chaotic motion often arises from a system with many degrees of freedom. 
A classical example is chaos in a fluid system (such as Rayleigh-Benard convection), where 
very high-dimensional chaotic motion should underlie at a molecular scale. A canonical 
answer for the condition to have low- dimensional chaos at a macroscopic level is given by 
separation of scales distinguishable from a microscopic level. Still it is not clear how such 
separation is possible, since chaos can lead to the amplification of a small-scale error. 

To address the question, we note that a certain coupled dynamical system shows 
some lower dimensional collective motion for macroscopic variables, while microscopic vari- 
ables keep high dimensional chaos. To characterize such collective motion, Lyapunov expo- 
nent at a macroscopic scale will be introduced, which specifies the growth rate of error at 
macroscopic variables. By studying the dependence of the exponent on the length scale and 
the system size, it is shown how the "collective chaos" is compatible with microscopic chaos, 
and how they are separated at the 'thermodynamic limit'. 

Here we adopt a 'heterogeneous' globally coupled map (GCM) with a distributed param- 
eter: 

e N 

x n+l (i) = (1 - e)fi{x n {i)) + j- T Yl fj( x n(j)) (1) 

JV 3=1 

where x n (i) is the variable of the z'th element (i = 1,2,3, •■ ■ ,N) at discrete time n, and 
fi(x(i)) is the internal dynamics for each element. For the dynamics we choose the logistic 
map fi(x) = 1 — a(i)x 2 , where the parameter a(i) for the nonlinearity is distributed between 
[a — ^p, a + ^r] as a(i) = a + Aa ^ JV - > . As a macroscopic variable, we adopt the mean 
field, 

1 N 

hn = ^J2fiMi)). (2) 

iV i=l 

in which the collective motion is contained. On the other hand, chaos of x n (i) is referred to 
'microscopic' here. 

Conventional GCM with identical parameters is given by Aa = 0, whose study has 
revealed clustering, chaotic itinerancy, partial ordering, and so forth 0. In particular, 
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study of collective dynamics has gathered much attention |]-|J . When the coupling e is small 
enough, oscillation of each element is mutually desynchronized, and the effective degrees of 
freedom increase in proportion to the number of elements N. Still, a macroscopic variable is 
found to show some kind of ordered motion distinguishable from noise, ranging from torus 
to high- dimensional chaos [Q-Q]. 

For instance, Fig.|l] gives a return map of the mean field dynamics of our model ([!]), which 
shows some pattern that may suggest low-dimensional chaos. Torus motion is also found 
by changing the parameters [|J. Here we try to characterize such collective dynamics, and 
search for low-dimensional collective chaos. 

First note that the conventional Lyapunov exponents are not relevant to the character- 
ization of collective motion. Indeed, for small e, all the N Lyapunov exponents of (1) are 
positive (whose values are close to the exponent of a single logistic map x — > fi(x)). No ex- 
ponent corresponding to the mean field motion is observed in this spectrum of N exponents. 
This apparent paradox can be resolved by noting the order of limit to define the Lyapunov 
exponent. In the calculation of the Lyapunov spectrum we take a limit of disturbance 
applied to the orbit. As long as we choose this limit first and then the thermodynamic 
limit (AT — > oo), the N Lyapunov exponents cannot characterize the collective motion ||. It 
is necessary to take the thermodynamic limit first and then the limit of disturbance scale, 
to characterize the collective dynamics. 

Since we are concerned with a system of large but finite N, the above order of limit implies 
that we have to keep the disturbance amplitude finite. To study such orbital instability, the 
finite-size Lyapunov exponent introduced by Vulpiani and his coworkers || is useful. It is 
given by 



where r is the maximum time such that \x' n — x n \ < A for trajectories x n and x' n starting 
from xq and x' Q = x$ + 5q respectively, while (•) is an average over the trajectories starting 
from different initial values. The length scale A can be considered as the scale of observation. 




(3) 
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We have measured the finite-size Lyapunov exponent for the macroscopic variable h n of 
our GCM ([!]). Here we perturb the orbit to give rise to a change from ho to h' = h Q + S (see 
the caption of Fig.^Jfor detailed description). In Fig.|2|, the finite-size Lyapunov exponent is 
plotted with the change of A. 

As long as the system size is finite, this finite-size Lyapunov exponent reflects not only the 
macroscopic motion but also the microscopic chaos. Indeed the exponent in Fig.[| changes 
with the scale A and no clear plateau (except for So — > 0) is visible. On the other hand, 
if low-dimensional macroscopic dynamics has a characteristic time scale separated from the 
microscopic dynamics, it will be possible to extract the growth rate of perturbation in the 
collective motion from the finite-size Lyapunov exponent for the macroscopic variable. To 
do so, we postulate the following assumptions that are expected to hold if the collective 
dynamics is low-dimensional chaos or on a torus. 

First note that in the limit, A — > and So — > 0, the finite-size Lyapunov exponent A,5 (A) 
converges to maximum Lyapunov exponent X m , which is determined by the conventional 
Lyapunov exponents for the microscopic variables x n (i) directly. 

Considering that the collective dynamics appears by coarse-grained macroscopic vari- 
ables, we postulate that there are length scales (in the phase space) A G [A m ,A c ], where 
the macroscopic variable is characterized by "collective Lyapunov exponent" Ac- Below 
A < A m the microscopic chaos dominates, while the orbit is out of the attractor (at a 
macroscopic level) for A > A^. To have low-dimensional collective dynamics, it is pos- 
tulated that Xc is independent of N (as long as it is large enough), and that A m should 
approach zero with iV — > oo while Ac remains finite. 

Based on the above assumptions, we can derive an approximate form of the finite-size 
Lyapunov exponent against the scale A. Let S n denote the distance from the original tra- 
jectory at time step n. For the scale A < A m , S n increases proportionally with exp (X m n). 
Hence t(A) = j- log — follows, independently of the collective dynamics. 

On the other hand, for the scale A m < A < A c , S n is given as S n oc exp (Xcn) for a chaotic 
case with Ac > 0, or S n oc n K for a torus case with k as a certain constant. Corresponding 



to each collective motion, r(A) and the finite-size Lyapunov exponent A,5 (A) are given by 



r(A) 



is log -t + i lo § If ( chaotic case ) (4) 

(£-)- + £]og/£-l (tornscase)' 



and 



A m A c log £ 



A m , : ,° 5 - (chaotic case) 



A C log^+A m log-£- 



A 5o (A) = <! " "1 og A — ^ . (5) 



*0 



(torus case) 



Fig.^ shows the behavior of Aa (A) for our GCM, with a fitting curve by Eq.(5). Here 
the parameters Ac and A m are obtained to fit the data for several values of 5q and A. The 
estimated exponent Ac is positive, and is much smaller than the Lyapunov exponent A m 
for microscopic chaos. Thus there is a regime in which the mean field dynamics shows low 
dimensional chaotic dynamics, although there is no clear plateau corresponding to Ac- The 
value Ac is smaller than Aj (A). 

It is then necessary to study iV dependence of A m , in order to confirm the existence 
of the low-dimensional collective motion. For this, it is convenient to transform Eq.(5) to 
remove Sq dependence of the data. For it, we define t(A) as t(A) = r(A) + j- log<5 , which 
characterizes the time for amplification of error from a certain scale independent of 8q. From 
Eq.(4), we obtain 

{T—- - log A + [ -ri -J- ) log A m (chaotic case) 
Ac ^ VA m a c ; > _ (6) 

(A) K + i log Am - 1 ( torus case ) 

Thus, the N dependence of A m appears as a shift of constant in t — A plot, while Ac or k 
is given by a slope in a suitable plot. 

In Fig.||, t is plotted as a function of A. As is shown in Fig.f|(a), the slope of the semi-log 
plot is independent of N. The Lyapunov exponent Ac, characterizing the collective motion, 
is given by the inverse of the slope, and is estimated as 0.02. On the other hand, A m , given 
by the shift of the plots, decreases with N, while Ac does not show significant change. Thus 
the scale for the collective motion A m < A < Ac increases with N. In Fig.|5], iV dependence 



of A m is plotted, which gives A m ~ whose form is expected from the central limit 
theorem. Hence the emergence of low-dimensional collective chaos at the thermodynamic 
limit is confirmed. 

We have also applied the present algorithm to the case with a collective torus motion. 
Fig.f|(b), (t — A plot), shows that k, the inverse of the slope, is 0.5, independent of N. 
Indeed this exponent 1/2 is expected from the diffusion process of phase on the torus. The 
decrease of A m with iV is also plotted in Fig.||, which again shows the expected decrease of 
A rn ~ ^= . Hence the collective torus motion is demonstrated. 

In this letter, we have proposed an algorithm to characterize the collective (chaotic) 
motion, and applied to it to a GCM. We have introduced collective Lyapunov exponent, 
to characterize the growth rate of perturbation in the collective motion. The microscopic 
chaotic motion exists at a small scale of the macroscopic variable, but such scale A m is shown 
to decrease as \/yN. Hence, the macroscopic motion is separated from the microscopic 
motion and the emergence of low-dimensional collective motion with N — > oo is confirmed. 

Existence of low-dimensional collective chaos in the presence of microscopic chaos has 
often been suspected JTD|. Indeed for a GCM with homogeneous elements (i.e., with Aa = 0), 
such low- dimensional collective chaos has not been observed so far. In Fig.[|(c), we have also 
applied our algorithm to this case. The separation of scales is not clear and the data cannot 
be fitted with (6). The shift of the plot gets smaller with the increase of N. At leaset A rn 



does not decrease as 1/ \JN JTT 



Our t — A plot provides a tool to distinguish low-dimensional collective chaos from high- 



dimensional one. In the former case, the plot shifts as log (Jl/N) with N, while for the latter 
case such shift is not observed. This distinction generally holds, even if the approximation 



to get (6) may not be very good JT3[ 



Our estimation of macroscopic motion is realized for a system subjected to microscopic 
chaos. It is expected that our algorithm can be applied even if we do not know the equation 
of motion, since the method of || is based on Wolf's algorithm |T2] developed for the estimate 
of Lyapunov exponents from experimental data. Thus, we hope that our method developed 

6 



in this letter is applicable to data obtained from experiments. 
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FIGURES 

FIG. 1. An example of return map for chaotic collective motion, ao = 1.92, Aa = 0.088, 
e = 0.1, N = 10 . Points (h n ,h n +\) are plotted over 3 x 10 4 steps after transient are discarded. 

FIG. 2. Aft (A) is plotted for the model © with a = 1.92, Aa = 0.088, e = 0.1(0), a = 1.9, 
Aa = 0.05, e = 0.098(x), a = 1.9, Aa = 0.05, e = 0.11(A), a = 1.69755, Aa = 0.0, e = 0.008(a), 
with N = 10 7 . Initial perturbation amplitude do is fixed at 1.0 x 10~ 7 . For computation, dis- 
placement H'q = ho + 5q is created by perturbing the orbit as x' (i) = xo(i) + 5q x a, where a 
is a random number in [—1,1]. Each point is obtained by averaging over 100 samples. Specific 
choice of this perturbation scheme is irrelevant to our results, as long as the collective variable is 
perturbed. Adopting the algorithm to be presented, the collective motion is shown to be torus 
(A), low-dimensional chaos (0 and x), and high-dimensional chaos (□). 

FIG. 3. The finite-size Lyapunov exponent A,5 (A) is plotted as a function of A for several 
initial displacement <5o- The exponent is computed from the average over 100 samples starting 
from different initial condition. The curve from Eq.([|) fitted to the data is also indicated. From 
the fitting, A m = 0.00086 and Ac = 0.02 are obtained, while A m = 0.41 is directly obtained from 
Eq.fllJ). a = 1.9, Aa = 0.05, e = 0.098, and N = 10 6 . Note that although A 5o (A) for small A is 
slightly smaller than A m , it approaches A m with 5q — » 0. 

FIG. 4. The normalized time steps i(A) are plotted for N = 10 4 , 10 5 , 10 6 , and 10 7 , with 
the fitted curves(^). (a) chaotic case (with a semi-log plot), for ao = 1.9, Aa = 0.05, e = 0.098. 
(b) torus case (with a log-log plot), for ao = 1.9, Aa = 0.05, e = 0.11. The maximum Lyapunov 
exponent A m =0.41 (a), 0.39(b) are obtained directly from our model(|l]). The parameters obtained 
by a least square fitting algorithm give Ac = 0.02. (a), k = 0.5(b). (c)high- dimensional case, which 
does not obey Eq.(||), (with a semi-log plot), for ao = 1.6962, Aa = 0, e = 0.008. In this case, 
while the return map shows some structure, t for N = 10 6 and 10 7 are not separated any more. 



10 



FIG. 5. The microscopic length scales A m are plotted as a function of N for several parameters. 
A m is obtained from the fitting indicated in Fig.||. The parameters are ao = 1.92, Ao = 0.088, 
e = 0.1(0, chaos), a = 1.9, Aa = 0.05, e = 0.098(x, chaos), a = 1.9, Aa = 0.05, e = 0.11(A, 
torus), 
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Fig. 1 T. Shibata & K. Kaneko 
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Fig. 2 T. Shibata & K. Kaneko 




Fig.3 T. Shibata & K. Kaneko 




Fig.4(a) T. Shibata & K. Kaneko 
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Fig.4(b) T. Shibata & K. Kaneko 




Fig.4(c) T. Shibata & K. Kaneko 




Fig. 5 T. Shibata & K. Kaneko 
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